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Abstract. Modeling of interacting galaxies suffers from an extended pa- 
rameter space prohibiting traditional grid based search strategies. As an 
alternative approach a combination of a Genetic Algorithm (GA) with 
fast restricted N-body simulations can be applied. A typical fit takes 
i^ ' about 3-6 CPU-hours on a Pentiumll processor. Here we present a par- 

allel implementation of our GA which reduces the CPU-requirement of a 
parameter determination to a few minutes on 100 nodes of a CRAY T3E. 
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Introduction 



A major problem in modeling of encounters of galaxies is the extended parame- 
ter space which is composed by the orbital and the structural parameters of the 
interacting galaxies. Traditional grid-based fitting strategies suffer from very 
large CPU-requirements. E.g. for a restricted 7-dimensional parameter space 
(an encounter of a disc with a point mass) and a resolution of only 5 values per 
i-^ I dimension, one needs 78125 models, or about 26 years of integration time on a 

GRAPESaf special purpose computer for a "complete" grid. More systematic 



O I search strategies like gradient methods depend strongly on the initial conditions, 

^— > . which makes them prone to trapping in local optima. An efficient alternative 

~^ \ approach are evolutionary methods and especially genetic algorithms (Holland 

1975, Goldberg 1989, Charbonneau 1995). In combination with fast (but not 
self-consistent) restricted N-body-codes (Toomre &: Toomre 1972) they allow for 
^ ' an efficient search in parameter space which can be used for both, an automatic 

H I search of interaction parameters (provided sufficiently accurate data are avail- 

able) and/or a uniqueness test of a preferred parameter combination (Wahde 
1998, Theis 1999). 

The basic idea of GAs is to apply an evolutionary mechanism including 
'sexual' reproduction operating on a population which represents a group of 
different parameter sets. All members are characterized by their fitness which 
quantifies the correspondance between simulations and observations. In order 
to determine the 'parents' two individuals are selected according to their fitness. 
These parents represent two points in parameter space. The corresponding pa- 
rameter set is treated like a 'chromosome', i.e. it is subject to a cross-over and 
a mutation operation resulting in a new individual which is a member of the 
next generation. Such a breeding is repeated until the next generation has been 
formed. Finally, the whole process of sexual reproduction is repeated iteratively 
until the population confines one or several regions of sufficiently high fitness in 
parameter space. 

Here we present a parallel implementation of our genetic algorithm. 
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Figure 1. 



left: parallelization scheme; right: speed-up vs. number 
■° Note the better performance, if the workload per 



iVpE of processors. Note the better performance, if the workload per 
processor (A^pop/-^PE) is kept constant (A^pop is the population size of 
192 here). 



2. Parallelization of the genetic algorithm 

We parallelized our code on the level of the GA by applying a 'master-slave' 
technique using message passing (MPI). On the master processor all the repro- 
duction operations are performed. It sends the individual parameter sets to 
the slave processors which do the N-body simulations and determine the corre- 
sponding fitness values (Fig. 1 left). The speed-up, i.e. the ratio of CPU-times 
ti and tjv used by a calculation with one or N processors, respectively, reaches 
values between 30 and 60 for 100 processors (Fig. 1 right). The deviation from 
the optimal gain is mainly caused by the spread of CPU-times required for the 
individual simulations. It becomes more important with decreasing population 
size. By parallelization the whole GA-fit using 10000 models reduces to less 
than 3 CPU minutes on a CRAY T3E. Thus, an 'interactive' analysis becomes 
possible. 
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